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Abstract 

NiTi porous materials with unprecedented superelasticity and shape memory 
were manufactured by self-entangling, compacting and heat treating NiTi wires. 
The versatile processing route used here allows to produce entanglements of 
either superelastic or ferroelastic wires with tunable mesostructures. Three di¬ 
mensional (3D) X-ray microtomography shows that the entanglement mesostruc- 
ture is homogeneous and isotropic. The thermomechanical compressive behavior 
of the entanglements was studied using optical measurements of the local strain 
field. At all relative densities investigated here 25 — 40%), entanglements 
with superelastic wires exhibit remarkable macroscale superelasticity, even af¬ 
ter compressions up to 25%, large damping capacity, discrete memory effect 
and weak strain-rate and temperature dependencies. Entanglements with fer¬ 
roelastic wires resemble standard elastoplastic fibrous systems with pronounced 
residual strain after unloading. However, a full recovery is obtained by heating 
the samples, demonstrating a large shape memory effect at least up to 16% 
strain. 

Keywords: NiTi porous materials, Entangled Monofilament, Architected 


* Corresponding author 

Email address: Laurent.0rgeas@3sr-grenoble.fr (L. Orgeas) 


Preprint submitted to Acta Materialia 


May 12, 2015 



Materials 


1. Introduction 

Due to their biocompatibility and potentially interesting mechanical prop¬ 
erties, porous NiTi materials are promising architected media, in particular 
for biomedical applications such as bone implants mm- Up to now, most of 
these materials were processed by powder metallurgy technologies, e.g. self- 
propagating high-temperature synthesis 00 , spark plasma sintering [5], hot 
isostatic pressing [6] and conventional sintering mmm- Many difficulties arise 
with these processing routes, including the growth of a brittle oxide skin and a 
difficult control of (i) the porous mesostructure, (ii) the microstructure and (iii) 
the thermomechanical properties of the NiTi solid phase. Intense research was 
devoted to overcome these difficulties, largely improving the resulting mesostruc- 
tures and biocompatibility mm- However, the macroscale mechanical proper¬ 
ties of NiTi porous media, such as their superelasticity and shape memory, are 
still lower than in bulk NiTi alloys. Furthermore, porous alloys suffer from large 
residual strains induced by stress concentrations that arise at grain boundaries 
and necks and lead to local plastic deformations and ruptures m- 

In parallel, there is a growing research activity on the development of new 
types of metal porous media made of self-entangled monofilaments mmm 
mmmm- These fibrous materials are often called ” entangled metallic wire 
materials” or ” metal rubber” because their mechanical behavior reminds that 
of elastomers. Compared to NiTi foams, manufacturing is simpler, consisting in 
cold/hot working treatments to create a significant amount of self-entanglements 
of the wire. Up to now, various metal wires have been used: steel mmm, 
A1 [12], Ti [IS] US] H] and Ni [15] . To strengthen their mechanical response, 
sintering or glueing operations can be used djnBmanaini. In all cases 
however, the mechanical behavior still exhibits limited reversible strain, typi¬ 
cally 0.01 — 0.1 in compression. This is mainly due to the low yield stress of 
the elastoplastic wires, which deform plastically during the macroscopic defor- 
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mation of the entanglements. To limit this effect, Courtois et al. m entangled 
high yield stress perlitic steel wires, but the resulting architecture was highly 
heterogeneous due to the important spring back effect induced by the stiff wire 
and the entanglements had to be confined inside a solid container, thus severely 
limiting the use of perlitic steel wires. 

In this context, the main objective of the present work is to process a novel 
homogeneous porous NiTi material that exhibits important superelastic and 
shape memory effects, i.e. with a strain reversibility largely above those ob¬ 
served with other NiTi porous media, as well as other entangled monofilaments. 
To circumvent the difficulties mentioned above, we revisited and adapted the 
manufacturing route of single-wire materials to the case of a commercial biomed¬ 
ical grade NiTi alloy. Moreover, tuning the processing parameters allowed us to 
tailor the mesostructure, superelasticity behavior and shape memory properties 
of these new NiTi porous materials. In the following, the methods used to pre¬ 
pare and characterize these materials are presented (Section [ 2 ]), together with 
both the properties of the shaped wires (Section [3| and their mesostructures 
(Section |4|. The macroscale thermomechanical properties of the entanglements 
are then investigated and discussed (Section [ 5 ]). 

2. Experimental procedure 

2.1. Material and sample processing 

The route used to manufacture NiTi porous materials is inspired from the 
work of Tan et al. [12] for ductile and elastoplastic Al alloy wires. As illustrated 
in Fig. [l] there are three steps to process homogeneous specimens: (i) shap¬ 
ing a spring from a straight wire using plastic deformation, (ii) entangling the 
spring into a spring ball, and (iii) setting the final shape and consolidating the 
entanglement by close die compaction. 

A standard biomedical-grade Ni-50.8at.%Ti supplied by Forth Wayne Metals 
in cold worked state with a diameter d = 500 /am was chosen for its capability to 
display a superelastic behavior at room temperature and a shape memory effect 
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after suitable short and low-temperature heat treatments. The spring shape 
in step (i) was achieved by rolling the wire around a 3 mm diameter threaded 
cylindrical bar, with a pitch close to the wire diameter. A cylindrical steel close 
die mounted on tension-compression testing machine (MTS DY34, maximum 
force lOOkN) equipped with a furnace was used for step (iii). The final samples 
have a cylindrical shape with a height Ho = 27 mm and a diameter Do = 20 
mm. Samples of various initial relative density, or fiber content, </>, ranging from 
0.27 to 0.36, were produced by starting with straight wire of lengths between 
11.66 m and 15.55 m. 

In its initial cold worked state, the NiTi wire is elasto-plastic with a high 
yield stress (~ 1800 MPa), pronounced strain hardening and a rather small 
ultimate tensile strain 0.04). These features are not compatible with the low 
yield stress and high ductility necessary to process homogeneous self-entangled 
materials, as already pointed out in the Introduction mmm- To overcome 
these difficulties, the deformed wires were subjected to short optimised heat 
treatments dung at 623 K for 2 min (followed by water quenching) after both 
steps (i) and (iii). 

The temperature and duration of the heat treatments were also chosen to en¬ 
sure that the wire was superelastic at room temperature jH HU, a state called 
state 1 hereafter. Therefrom, the entanglement thermomechanical properties 
could be tailored using additional heat treatments. For example, the superelas¬ 
tic domain could be shifted to higher temperatures by heating the samples at 
723 K during 1 h. The resulting new state of the entangled wire, which is now 
ferroelastic at room temperature, will be called state 2. 

2.2. Mesostructure characterization 

Information about the entanglement mesostructures was obtained by scan¬ 
ning two typical samples, i.e. with (j)= 0.3 and 0.36, inside an X-ray laboratory 
microtomograph (RX solution, 3SR Lab, Grenoble, France). The scans were 
obtained with 1200 2D radiographs onto a 1914 x 1580 pixel flat panel (leading 
to a voxel size = 23 x 23 x 23 /rm 3 ), a scanning intensity of 200 /iA, a tension 
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of 150 kV, and an exposure time per radiograph of 0.25 s (to restrain the noise, 
an average of six radiographs per 2D image was used). To reduce beam harden¬ 
ing effects, a 1 mm Cu filter was positioned between the X-ray source and the 
sample. After reconstruction followed by filtering and thresholding operations, 
3D binarized images of the samples were obtained. An example is shown in Fig. 
|2|a). The images were used to analyse: 


• the local fiber content profiles along the height and radius of the samples 
(Fig. gc-d)). 

• the pore size distribution, by using morphological operations from a hexag¬ 
onal structural element with an eight voxel size (Fig. |2jf)) [23| [24] . 


• the anisotropy, connectivity and mean curvature of the wire, by follow¬ 
ing the method proposed by Latil et al. [25] and Orgeas et al. ‘26\. As 
illustrated in Fig. [2jb) , the centerline of the wire was first detected by 
(i) estimating the 3D Euclidian distance map of the fibrous phase, (ii) 
thresholding the resulting grey scale image at an intermediate grey value 
to eliminate contact zones from the image, (iii) skeletonizing the thinned 
wire, (iv) smoothing the resulting skeleton to obtain its centreline. The 
skeleton was discretized in elementary segments. The segment length was 
chosen equal to 2d, twice the fiber diameter. Each segment i has a curvilin¬ 
ear abscissa and unit tangent vector, respectively noted Si and p^. For the 
two scanned samples, i.e. with </>= 0.3 and 0.36, the number of segments 
N was 12960 and 15552, respectively. 


The mean orientation of the centerline was determined by computing the 
second order fiber orientation tensor a m (see Figge)): 

N 


a 4ep. 


N 


(1) 


i=1 


The variation of tangent vector along the curvilinear abscissa Si of the 
wire was used to compute the local mean curvature Ki (see Fig. [2|g)): 

11 d P 11 


ds 


( 2 ) 
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Self-contacts were defined as pairs of elementary segments with a distance 
less that the fiber diameter. Each contact j was ascribed an orientation 
unit vector n j normal to its contact plane. The mean contact orientation 
was then estimated from the second order contact orientation tensor B 


(see Fig. |2je)): 



(3) 


where M is the total number of elementary contacts. We should note 
that a single contact may be composed of several elementary contacts 
that link neighboring segments. We thus distinguish between elementary 
contacts as used in the above equation, whose number depends on the 
choice of discretization length, and full contacts that were obtained after 
agglomeration of elementary contacts between neighboring segments. In 
the following, we will note n, the number of full contacts per unit fiber 
length. Their length l c was also computed and analyzed (see Fig. |2jh)). 

2.3. Thermomechanical characterization 

The local stress-free forward and reverse martensitic transformations of the 
processed entanglements were analyzed by cutting small wire pieces inside sam¬ 
ples and subjecting them to DSC measurements. Straight wires that had under¬ 
gone the same thermal history as the entanglements were also tested. The tensile 
mechanical behavior was assessed at various temperatures ranging from 293 K 
to 353 K, subjecting the wires to load-unload cycles up to a maximum tensile 
strain of ~ 0.08 at a strain rate of 10 -3 s -1 (testing machine Gabo Eplexor 500 N 
equipped with a furnace). Furthermore, the processed entanglements were sub¬ 
jected to compression cycles at various temperatures and strain rates ranging 
respectively from 293 K to 353 K and from 4 10 -4 s -1 to 3.7 10 -2 s -1 . The sam¬ 
ples were placed inside a compression rheometer with parallel plates equipped 
with an inner heating system [28] and mounted on a mechanical machine (MTS 
DY26, maximum force 100 kN). During the tests, the recorded compression 
force, F, and the current sample height, F, were used to estimate the nominal 
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compression stress a = AF/ttDq and global axial compression Hencky strain 
e g = \il(H/Hq). In the following, the reference testing temperature and strain 
rate are set to 293 K and 4 10 -4 s _1 . 

To avoid errors on the strain measurements induced by edge effects such as 
friction between the sample and the plates or inhomogeneities on the outside 
of the samples, a local strain measurement was developed using images of the 
sample taken during the tests with a CCD camera (Jai Pulnix RM-4200GE, 
spatial resolution 2048x2048). As illustrated in Fig. [3j a segmentation of the 
grey scale images (Fig. [3ja,b)) together with a particle tracking procedure 
(developed with ImageJ and Matlab) were used to obtain the coordinates of 
the centers of mass of the enlightened parts of the fiber (Fig. Hc,d)). By 
following the image sequences, we defined a gauge zone in the center of the 
samples, within which the strain could be assessed. The enlightened parts of the 
fibers that were as close as possible from the borders of this zone (in the initial 
configuration) were automatically detected and followed during the compression 
cycles to determine the initial (ho) and current ( h ) heights of the gauge region, 
from which we determined the local Hencky axial strain e = \n(h/ho). 

3. Wire properties after processing 

As shown in Fig. |4ja), within the investigated temperature range, DSC mea¬ 
surements for the entangled wires in state 1 exhibit a single transformation peak 
upon cooling and heating, presumably associated with forward A-R (Austenite 
to R-Phase) and reverse R-A martensitic transformations These peaks are very 
flat, with latent heats A H A ~ R ~ A H R ~ A ^ 2 J g -1 . Conversely, in state 2 , 
the wires exhibit two stage A-R-M (Austenite to R-phase to Martensite) and 
M-R-A transformations with pronounced transformation peaks and latent heats 
(A h a ~ r ~ m ~ AH M ~R-A ^ 16 J g -1 ). Upon heating, the transformation end 
temperature, Af, is close to 328 K. 

As expected from the DSC measurements, the room temperature wire ten¬ 
sile behavior in state 1 is superelastic (see Fig. 0b)): upon loading, after an 
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initial apparent elastic regime, the stress-strain curve shows a high R-M trans¬ 
formation stress, cr tr « 550 MPa, and a transformation plateau within a strain 
range A e tr ~ 0.05. Unloading leads to the reversed process, with a pronounced 
stress hysteresis and a small residual strain. Conversely, wires in state 2 are fer- 
roelastic, with a large residual strain upon unloading due to the unrecoverable 
strain produced during loading by the stress-induced A-R-M and A-M forward 
transformations. Note however that the transformation strain may be recovered 
by stress-free heating above Af. Also, as seen in Fig. Qc), the transformation 
stresses recorded in both states at the onset of the stress-induced forward trans¬ 
formation, cr tr , follows a Clausius-Clapeyron like relation, with a slope close 
to 5.6 MPa.K -1 , in agreement with values commonly reported for similar heat 
treated NiTi alloys [29l I2T1 [22 . 

4. Mesostructure of the entanglements 

Fig. [^summarizes the information gathered on the entanglement mesostruc¬ 
ture. We focus here on 2 representative samples with mean relative densities 
0 =0.3 and 0.36. Fig. |2|a) shows a cut through a tomographic reconstruction, 
where we see a complex architecture that may be described at two levels: at 
short range, the structure is made of a quasi-ordered spring, but at longer range, 
the spring is entangled with itself in a disordered fashion. As illustrated from 
the two profiles of the normalized relative density in Fig. §c-d) (0* = 

< j) and 0 being the local and mean relative densities), the fiber content inside 
the samples is close to homogeneous, except slight decreases of 0* close to the 
surface of samples. This is in contrast with systems obtained after entangling a 
straight wire where the density strongly increases on the outside of the samples 
p~4] : in this situation, the reason is that the fiber decreases its curvature by 
moving to the outside. Here, this relaxation is not possible because the high 
tortuosity of the spring leads to the formation of numerous contacts in the spring 
ball (step (ii) in Fig. [l]) that efficiently lock the fiber. 

The pore size distribution shown in Fig. §f) exhibits a broad peak at 


~ 1300 /i m for </> = 0.3. As expected, increasing the fiber content to 0 = 0.36 
creates smaller pores, with a sharper peak shifted down to « 750 /im. The peak 
pore size may thus be tuned by simply changing the length of the initial wire 
(and thus the fiber content 0) without changing any of the other processing 
parameters. This will also affect the mechanics of the entanglements (see next 
Section). It is also interesting to notice that the pore size distribution could 
also be shifted down at fixed fiber content (and presumably fixed mechanical 
properties) simply by reducing both the wire diameter and the diameter of the 
threaded bar used to coil the springs. Thus, combined with the open-pore nature 
of the entanglement, the resulting mesostructure would fit values commonly 
admitted for bone implants, i.e. mean pore size between 100 and 600 fi m [l|. 

Distributions of the mean curvature are reported in Fig. [2|g). They display 
narrow peaks near 0.55 mm -1 and 0.52 mm -1 for </> =0.3 and 0.36, respec¬ 
tively. These values are close to the principal curvature Rq ~ 0.57 mm -1 of the 
non-entangled spring. Thus, the moderate stretching and folding of the spring 
induced during steps (ii) and (iii) of the manufacturing process do not induce 
noticeable changes of this mesostructure parameter. 

The diagonal values of the orientation tensors A and B are nearly identical 
and close to 1/3 (see the data given in Fig. [2^e) for the sample with 0 =0.36), 
confirming that on average, the entanglement is isotropic, without preferred 
orientation for the wire centerline or its self-contacts. 

The estimated number of full contacts per unit fiber length, n, is close to 
0.40 mm -1 and 0.50 mm -1 for </> =0.27 and 0.36, respectively. In other words, 
the two scanned samples exhibit respectively 2608 and 3874 full contacts. This 
emphasizes both the very high connectivity of the entanglements and its increase 
with fiber content 0. Besides, the measured ratio r = n(0 = 0.3 )/n(<fi = 0.36) is 
close to 0.83. It is worth mentioning that this ratio is consistent with the ratio 
r predicted by the excluded volume theory for isotropic networks of straight 
fibers with circular cross section mi eh mam]: 

r(<t> = <t>i/<t>2, (4) 
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which yields r = 0.81 for =0.3 and (j )2 =0.36, again emphasizing the effect of 
the fiber content on the entanglement connectivity. 

Lastly, as shown in Fig. the distributions of contact lengths are wide 

with a predominant contact population marked by a distinct peak, but also a 
distribution of short contacts, down to zero length. It is also worth noting that 
the predominant contact length respectively peaked close to 400 fi m and 700 /im 
for 4> =0.3 and 0.36: increasing the fiber content yields to an increase of both 
the number of full contacts (see previous point) and their length. 

5. Thermomechanical behavior 

5.1. Typical results for samples in state 1 

Fig.0 shows typical stress-strain curves obtained during compression load- 
unload cycles at the reference testing temperature and strain rate. Here the 
stress is shown as a function of both the global and local strains. At the sample 
scale, Fig. [3] shows that the deformation proceeds without pronounced friction 
effects at the compression plates since the initial cylindrical shape is preserved 
even at finite deformation. However, we have seen in Fig. [2] that the relative 
density is not fully homogeneous, in particular near the top and bottom of the 
samples. We may thus expect more pronounced deformations in this regions. 
Indeed, we find in Fig. [5] that the global strain measured at the level of the 
structure is noticeably larger than the local strain, which is measured in a 
gauge region near the center of the sample. We ascribe this difference to larger 
deformations near the top and bottom of the samples. For this reason, we will 
consider only the local strain in the following. 

As shown in Fig. [5j the loading curve is practically recovered during succes¬ 
sive cycles. The fibrous entanglement thus exhibits no stress accommodation, 
which is rarely observed in non-bonded fibrous materials, where irreversible 
rearrangements of the fiber segments, or, in case of elastoplastic (resp. brit¬ 
tle) fibers, plastic deformation (resp. rupture) of fibers, occur during cycling 
021 GHEE]. This marked stability of the entangled structure arises from the 
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numerous self-contacts and the fact that during deformation, sufficient internal 
back stresses arises in the wire to reverse during unloading the rearrangements 
induced upon loading, at least at the macroscopic level of the structure. 

The stress-strain loading curve exhibits a marked strain hardening, even at 
elevated compression strains, which is interesting for structural applications. 
This was observed independently of the state, fiber content, temperature and 
strain rate. This phenomenological feature is different from foams [35] . which 
exhibit stress plateaux. The reason is the continuous formation of new contacts 
during compression, which stiffen the structure. The stress-strain curves in Fig. 
[5] start with a low slope, at least during the first cycle, which is typical of non- 
bonded fibrous systems where the fibers rearrange to adapt to the applied strain. 
Then, the strain-stiffening rate increases rapidly up to about 5 = 0.1. This be¬ 
havior is related to a consolidation of the structure mum urn cm, the principal 
underlying mechanisms being a rearrangement of the fibrous mesostructure, an 
increase of the number and/or surface of fiber-fiber contacts, which in turn re¬ 
strains the motion, rearrangement and deformation of fiber segments [39, 155] . 
For strains above about 0.1, the stress increases less rapidly, which is unusual for 
compressed fibrous materials made of elastic fibers. In the case of elastoplastic 
(or brittle) fibers, this phenomenon is related to the plasticity (resp. rupture) of 
fibers during loading Banana- For NiTi wires, it should rather be ascribed to 
the forward martensitic transformation induced in the vicinity of contact zones 
where the wire is significantly bent. 

Fig. 0 also shows upon unloading a nearly full recovery of the imposed 
strain, even when the maximum local strain reaches 0.3 (corresponding to a 
global strain « 0.4). The strain recovery is accompanied with a noticeable hys¬ 
teresis. To the best of our knowledge, such macroscale superelastic behavior 
has never been observed at elevated strains in porous/fibrous metals, even with 
similar mesostructures. Indeed, Tan et al. [12] and Zhang et al. 16 observed 
a nearly full compression strain recovery for entangled elastoplastic Al alloy 
monofilaments but with a maximum applied global strain of 0.13 only to avoid 
wire plasticity. The present macroscale superelastic behavior is induced by the 
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elastic energy stored in the wire upon loading and released upon unloading. It 
is also probably related to stress-induced reversible martensitic transformations 
taking place during the unloading phase. Both phenomena offer sufficient driv¬ 
ing force to the wire to reverse the deformation mechanisms that occurred upon 
loading, despite significant mesostructure rearrangements and friction forces at 
contacts. 

5.2. Influence of the strain rate and testing temperature 

Fig. © shows the evolution of 00 . 2 , the stress at e =0.2, as a function of 
the strain rate (Fig. [6ja)) and testing temperature (Fig. j6^b)). The strain 
rate has a weak influence, with an increase of a 0.2 by less than 20% when the 
strain rate is increased by 3 orders of magnitude. This weak increase may be 
related to the occurrence of thermomechanical couplings during stress-induced 
transformation in the wire mmm- Fig. [6^b) shows the evolution of 00.2 for 
a sample compressed at e = 4 10 -4 s -1 at temperatures from 293 to 353 K. 
Obviously, the effect of the temperature is also limited, with a slope of about 
0.05 MPa.K -1 , more than two orders of magnitude smaller that for a single wire 
(see Fig. [4|. Thus, results emphasized in Fig. [6ja) and (b) lead us to conclude 
that during the investigated compression tests, a very small volume fraction of 
the wire is transformed. 

5.3. Influence of the fiber content and the maximum compression strain 

Fig- 0 shows the influence of the initial fiber content and of the maximum 
compression strain imposed in a given load-unload cycle, on the mechanical 
behavior of the entanglements. As illustrated in Fig. [7](a) , the higher the fiber 
content, the higher the stress level upon loading, the larger the stress hysteresis 
and the larger the tangent moduli of stress-strain curves. These trends are often 
observed in non-bonded fibrous materials [36j [34] 37J [14] • 

The evolution of stress levels with fiber content is more precisely emphasized 
in Fig 0b). For the present range of fiber contents, the increase of 0 - 0.2 is 
limited: from </> = 0.27 to 0.36, 0 - 0.2 raises only from 3.6 to 6 MPa. This trend 
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is directly related to the increase of the number and length of self-contacts, as 
emphasized in the last Section (Fig. [ 2 ]): the larger the number and length of self¬ 
contacts, the more difficult the rearrangements and the higher the macroscale 
stress [39j [36] EH |25j [38]. 

The maximum hysteresis stress cr™y S x of a given load-unload cycle was esti¬ 
mated by measuring the maximum stress gap between loading and unloading, 
as illustrated in Fig. [7j^c). As shown in Fig. [7|d), <J™y S x is approximately twice 
larger when </> =0.36 than when </> =0.27. Besides, a nearly linear increase of 
v™™ with s 1710 * 00 is also noticed. These two trends can be explained by three 
underlying phenomena. First, the hysteresis associated with forward and re¬ 
verse stress-induced martensitic transformations increases with the maximum 
strain applied during a cycle [42] 03[ 04]. However, this phenomenon should 
be of second order given the small volume fractions of transformed wire es¬ 
timated above. The second possible mechanism is the rearrangement of the 
fibrous mesostructure: although the entanglements exhibit a superelastic be¬ 
havior without marked residual strain at the macroscopic level, the rearrange¬ 
ments at the wire level follow different paths during loading and unloading (see 
next Subsection). Lastly, there is the friction induced by forces arising at self¬ 
contacts, which is closely linked with the previous phenomenon. 

Tangent moduli just before, dcr + /d£, and just after, dcr _ /d6, the return point 
of a given load-unload cycle (see Fig. §b)) were estimated (dcr /&£ estimates 
the Young modulus E of the entanglements, as often proposed for foams or 
fibrous materials). Their evolution with £ ma:r are reported in Fig. He). As 
expected, da + /d£ is much lower than dcr - /d e, showing that, as expected, non¬ 
elastic mechanisms, such as wire rearrangements or friction as mentioned above, 
are involved during loading. Moreover, the measured moduli increase with s max , 
which could be correlated to a possible increase of the actual fiber content during 
loading and/or the increase of the number/surface of fiber contacts (induced by 
the possible increase of the actual fiber content or by the rearrangement of the 
entanglements). Furthermore, a noticeable increase of E with the initial fiber 
content </> is noticed: this is directly connected to the increase of the number of 
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fiber contacts with </> (see section [ 4 J [36]). 

Lastly, the Young moduli estimated from da~/ds in Fig^b) strongly de¬ 
pend on the fiber content and less on the applied strain. We find Young mod¬ 
uli broadly between 100 and 350 MPa. These are lower than often reported 
for metal foams with similar relative densities (in general above 1 GPa [35] 5. 
Nonetheless, the Young moduli found here are within the range of 10-3000 MPa 
estimated for trabecular bones nans ns, and very close to values given for 
porous NiTi materials used for biomedical implants m- The present architec¬ 
ture may therefore be a potential candidate for biomedical applications, with 
the very interesting advantage of being superelastic and highly adaptable to 
applied loads. 

5.^. Hysteresis and Discrete Memory 

The hysteretic behavior of the entanglements is further explored in Fig. [8] 
This stress-strain diagram was obtained using a specifically designed cycling. A 
first major load-unload loop noted 0-5-0 in Fig. [8] was performed from e = 0 
to e max = 0.3 then back to e = 0. Then, a partial loop cycling was started at 
the return point 1 and continued following the numerical order of the return 
points marked in Fig. [8] i.e. 1-2-3-4-5. We have also reported in Fig. [8] images 
extracted from the video recorded during the experiment. In Fig. [8^d), we 
computed the difference 3-3” from the image of the sample at the return points 
3 and 3” (Fig. (8^c)), i.e. two points with the same strain 6 but a different 
mechanical history. In this image, the larger the difference at a given pixel, the 
brighter the pixel. The same procedure was carried out between points 3 and 
3’ (Fig. [8|f)) and 1 and V (Fig. [8jg)). The grey scale histograms deduced from 
these differences are shown in Fig. |8jb). 

As clearly evidenced from Fig. [8ja), the stress-strain curve of the section 4-5 
passed through all the previous return points, thus closing all the incomplete 
partial loops, 3-4-3 and 1-2-1 initiated during the pathway 1-2-3-4. Fig. j8^d) also 
reveals some mesotructure differences between points 3 and 3”, although these 
two configurations have the same compression local strain. Such a difference is 
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also emphasized with the rather flat grey scale histogram of Fig. |8jb) . Both 
figures highlight the rearrangement of the fiber segments between the loading 
and unloading phases of a cycle. By contrast, images 3 and 3’, as well as 1 and 
V nearly coincide, despite the complex pathways used to go from 3 to 3’ and 
from 1 to I’. 

Similar behavior was observed in superelastic or ferroelastic bulk NiTi or 
CuZnAl alloys [43l 1 44] . As for these systems, the present entanglements fol¬ 
low the concept of “erasable micromemory” also known as the “return point 
memory” or “discrete memory” effect @3 031021. The re-assembling of the 
return points along the path 4-5 suggests that the mechanical behavior of the 
entanglements is ruled by a part of its previous deformation history. More pre¬ 
cisely, once a subloop is initiated, the entanglement ” remembers” the subloop 
starting/return point, namely its reference stress-strain state and it also remem¬ 
bers its reference mesostructure until the subloop is closed. Moreover, once the 
subloop is closed, it is “forgotten” in the way that it does not affect the par¬ 
ent loop within which it was performed. This is illustrated along the section 
4-5 after closing subloop 3-4-3 at return point 3: the stress-strain curve follows 
its parent loop 1-2-1 as if no subloop had been carried out before. Similarly, 
the stress-strain curve follows the major loop 0-5-0 as if subloop 1-2-1 and 3- 
4-3 had been erased from the entanglement mechanical history. This behavior 
could be potentially interesting for practical applications. From a fundamental 
standpoint, it raises questions about the mechanics of such non-bonded entan¬ 
glements that are able to memorize some of their mesostructure configurations, 
even after complex trajectories and rearrangements. 

5.5. Hysteresis and Energy Dissipation 

The observed hysteresis is related to dissipative mechanisms. To estimate 
these effects, the energy loss Eh ys induced during (sub)loops was computed and 
compared with the macroscale elastic energy E e i stored and released during such 
sequences. For each closed (sub)loop, Eh ys was estimated from the following 
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expression: 


Ehys — 



(5) 


The elastic energy E e i was calculated by assuming that the overall stress pro¬ 
duced by the entanglement can be seen as the sum of two contributions: a hys- 
teretic stress <Jh ys and a hyperelastic stress (J hypi defined as the middle (neutral) 
line of the major hysteresis loop (see Fig. [TJb) for an illustration) [47j 51, 16]. 
Hence, for a (sub)loop carried out between two strain states E\ and e 2 , E e \ is 
calculated as: 

pS 2 

E e l = / Ghyp^-S. ( 6 ) 

J £1 


The loss factor r/, i.e. 

m : 


the energy dissipation coefficient, is then computed as 


T] = 


Ehys 
nE e i ' 


(7) 


Fig. § a) reports the evolution of p with the strain magnitude A 5 = £2 — £1 
deduced from stress-strain diagrams such as shown in Fig. [ 5 ] {s\ =0 and 
£2 = £ max ). We can see that 77 depends on 0 than on £ max . As illustrated by 
the material map plotted in Fig. 0 b) , the values of 77 found here between 0.25 
and 0.3 are higher than most of metallic foams, because of the increased friction 
at contacts as well as the forward and reverse stress-induced martensitic trans¬ 
formations. Besides, with a much higher Young modulus, the entanglements 
have loss factors as high as the best polymeric flexible foams. Finally, the pro¬ 
cessed materials display loss factors similar to other entangled mesostructures 
mm whereas they can be deformed reversibly at higher compression strains. 
Such architected structures are thus potentially relevant design solutions for 
load bearing and damping components. 


5.6. Samples in state 2: ferroelasticity and shape memory effect 

By simply changing the heat treatments during sample processing and going 
from state 1 to state 2 , it is possible to tune the thermomechanical behavior 
of the entanglements, as illustrated in Fig. [lOj The two stress-strain curves 
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in Fig. 10 'a) were obtained with samples in state 2 compressed at two differ¬ 
ent maximum strains at the reference temperature and strain rate (</> = 0.3). 
The mechanical response is drastically modified compared with the stress-strain 
curves in state 1 (see Fig. [7](a) ). First, the stress levels are much lower. Second, 
the superelastic behavior of state 1 is replaced by a ferroelastic behavior with 
large residual strains. The residual strain reaches 0.08 and 0.16 for maximum 
imposed strains of 0.2 and 0.27, respectively. From a phenomenological stand¬ 
point, the shape of stress-strain curves is close to those observed for similar 
entanglements with elastoplastic wires m ■ However, for the present NiTi wire, 
the shape of the stress-strain curve is not controlled by the wire ductility but 
rather by the low value of the transformation stress required to induce marten¬ 
site within the wires (see Fig. |4](b)). Thus, in state 2 , the energy stored in the 
wire upon loading is not high enough to overcome the configurational and fric¬ 
tion forces required to push back the rearrangements and the wire deformation 
is highly irreversible upon unloading. 

The large residual strains induced at room temperature can be fully recov¬ 
ered simply by a stress-free heating up to 353 K, as evidenced in Fig. [lQ](b) : the 
samples come back to their initial shapes upon heating, thus exhibiting a large, 
full and thus unprecedented shape memory effect. Upon heating, the energy 
stored in the transformed zones of the wire is progressively increased according 
to the Clausius-Clapeyron relation and become sufficient to reverse the trans¬ 
formation in these zones, pushing back the rearrangements induced during the 
deformation at room temperature. 


6. Concluding Remarks 

By entangling a single NiTi spring followed by a compaction in a close die, we 
were able to manufacture homogeneous porous NiTi materials with controlled 
mesostructures and physical properties. The manufacturing method is versatile 
enough to tune easily the process-induced mesostructures. For examples, by 
changing the compaction rate or the initial wire length, samples with different 
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fiber contents may be obtained. Also, a simple scaling of the spring dimensions 
would allow to tailor the pore size distribution and changing the shape of the pre¬ 
compacted entanglement may also lead to orientated specimen with anisotropic 
properties. Also, compared with other methods used to produce NiTi porous 
materials, another advantage of the present processing route is the use of only 
low temperature heat treatments. This avoids the formation of an oxide skin, 
which is detrimental to the mechanical properties of NiTi. Another advantage of 
the present method is that the temperature and duration of the heat treatments 
can be tuned to optimize the thermo mechanical properties of the wire. An 
illustration was given here, where, by using higher temperature treatments, we 
switched from a superelastic to a ferroelasticity behavior with shape memory 
effect. 

The thermomechanical properties of the entanglements were analyzed based 
on compression experiments. Although other types of mechanical loadings, 
such as tension or shear, should be considered to complete the experimental 
database, the results obtained here are interesting and promising. Entangle¬ 
ments processed in state 1 exhibit a remarkable superelastic behavior with 
negligible residual strains even after repeated cycles. This behavior is weakly 
strain-rate and temperature dependent. The nearly full reversion of the macro¬ 
scopic strain observed upon unloading is explained by the unique association of 
superelastic wires together with their pronounced entanglement without strong 
physical bonds, allowing a reversible deformation and rearrangement of the wire 
during the compression without plastic deformation or rupture. This was not 
observed with other metal foams, nor with other entangled metallic monofila¬ 
ments. Combined with the relatively high stress levels and stiffnesses observed 
during compression (both of them could be tailored by increasing the fiber con¬ 
tent), this material could be used as stiffer substitutes to standard elastomeric 
foams. The Young modulus also corresponds to values obtained for trabecular 
bones, so that combined with their porous mesostructure, these entanglements 
could also be relevant solutions for bone implants or scaffolds. 

The superelastic behavior is accompanied with a stress hysteresis, which in- 
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creases with imposed strain and fiber content. Such a hysteresis is due to strain- 
induced rearrangements of the entanglements, to their associated friction forces, 
and to stress-induced forward and reverse martensitic transformations occurring 
in the wire. Likewise, the observed hysteretic behavior exhibits an interesting 
discrete memory, which needs to be clarified at the mesoscale and should also 
be taken into account in macroscale constitutive equations dedicated to model 
the mechanics of such materials. Lastly, we have shown that this pronounced 
hysteresis provides the entanglements with good damping properties, so that 
combined with their rather high stiffness, these entanglements can be viewed as 
suitable solutions for load bearing and damping structural components. 

Switching from a superelastic to a ferroelastic state profoundly changes the 
mechanical behavior of the entanglements, leading to a marked degradation of 
both the stress-levels and strain recovery. This is directly related to an appar¬ 
ent plastic deformation of the wire, induced both by martensitic transformation 
and martensite reorientation. We have shown that the corresponding macroscale 
residual strain, even at levels as high as 0.16, can surprisingly be fully recovered 
upon stress-free heating. This full shape memory effect has never been observed 
with other NiTi porous materials. Once again, it is related to the unique as¬ 
sociation of NiTi wire with a high level of entanglement free of strong physical 
bonds. 

To optimise the architecture of these materials and to better understand 
both their remarkable macroscale properties and their underlying complex de¬ 
formation meso-mechanisms, further and finer investigations are needed. In 
particular, performing similar experiments with 3D in situ observations of the 
entanglement mesostructures, e.g. by using X-ray microtomography mmm 
and/or using discrete/finite element simulations at the fiber scale [52l l53l 154] 
would undoubtedly bring relevant elements towards this end. This work is in 
progress. 
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Figure 1: Illustration of the various shape setting operations followed to produce a homoge¬ 
neous NiTi entangled single-wire material. 
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Figure 2: 3D view (a), centerline (b) and second order fiber and contact orientation tensors 
A and B (e) of a NiTi entanglement with a mean fiber content </> = 0.36, together with 
fiber content profiles <j>* = 4>/4> along its radius (c) and height (d). The pore size (f), mean 
curvatures (g) and contact lengths (h) distributions are also shown for both scanned samples. 
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Figure 3: Illustration of the particle tracking method used to assess the local compressive 
strain: grey scale (a,b) and segmented (c,d) pictures of the sample before (a,c) and during 
(b,d) compression. 
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Figure 4: DSC curves of small portions of the entanglements in states 1 and 2 (a), tensile 
stress-strain curves of two straight wires subjected to same thermal history as the entan¬ 
glements (b), and evolution of the forward stress-induced transformation stresses cr tr with 
temperature (c). 
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Figure 5: Typical stress-strain curves obtained during compression load-unload cycles for a 
sample a sample in state 1 (0 = 0.36, reference testing temperature and strain rate). 


(a) (b) 



Figure 6: Evolution of the compression stress cro.2 with strain rate e (0=0.3 and T=303K) 
(a) and testing temperature T (e=4 10 -4 s _1 ) (b). 
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Figure 7: Influence of the fiber content </> and maximum compression strain £ rnax on stress- 
stress curves (a), stress ctq .2 (b), maximum stress hysteresis (c-d) and tangent moduli 

dcr + /d£ and dcr _ /de (e). 


29 















(f) 


(g) 

( , 4 

i'\ ; . / ‘j/. . \ 


.i v ‘ ■ 

i, •, ■ ^ \ _ 


'■ - . ' d 

!? ■■ ■ 


; ^ 




Vi . v. 


> . ' .a ■-/ - 


3 - 3 ' 1 - 1 ' 


Figure 8: Complex cycling 0-1-2-3-4-5 performed on an entanglement to illustrate its discrete 
memory both on stress-strain diagram (a) and on the mesostructure (b-g): images taken at 
points 1, 1’, 3 (c), 3’ and 3” (e) were used to compute the image differences 3-3” (d), 3-3’ (f) 
and 1-1’ (g), the grey scale histograms of which are plotted in (b). 
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Figure 9: Evolution of the loss factor 77 with strain magnitude £ rnax for two fiber contents 
4> (a). Material map of the loss factor 77 versus the Young modulus (from CES [50))) (b): 
metallic entanglements are represented by rectangles, with the present NiTi architectures in 
black. 




Figure 10: Stress-strain curves (a) and strain-temperature curves (b) showing the ferroelastic 
behavior of the entanglements in state 2 and its associated shape memory effect observed upon 
stress-free heating after the load-unload sequence shown in (a) (0=0.3). 
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